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Abstract 

We give an asymptotic expansion {the higher Stirling formula) and an 
infinite product representation {the Weierstrass product formula) of the 
Vigneras multiple gamma function by considering the classical limit of the 
multiple g-gamma function. 



1 Introduction 

The multiple gamma function was introduced by Barnes. It is defined to be 
an infinite product regularized by the multiple Hurwitz zeta- functions [Q, Q, 
[Q, Q. After his discovery, many mathematicians have studied this function: 
Hardy Q] , ||^ studied this function from his viewpoint of the theory of elliptic 
functions, and Shintani applied it to the study on the Kronecker limit 

formula for zeta- functions attached to algebraic fields. 

In the end of 70 's, Vigneras ||2^ redefined a multiple gamma function to 
be a function satisfying the generalized Bohr-Morellup theorem, Furthermore, 
Vigneras , Voros and Kurokawa jl^ , |^ , , (l5| showed that it plays 
an essential role to express the Selberg zeta-function and the determinant of 
Laplacians. 

As we can see from these studies, the multiple gamma function is a funda- 
mental function for the analytic number theory: See also However we 
do not think that the theory of the multiple gamma functions has been fully 
explored. 

On the other hand, the second author of this article introduced a g-analogue 
of the Vigneras multiple gamma functions and showed it to be characterized by 
a (j-analogue of the generalized Bohr-Morellup theorem . 

*To appear in the proceedings of the workshop " Quantum groups and Quantum spaces" 
at the Banach center 
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In this article, we will establish an asymptotic expansion formula (the higher 
Stirling formula) and an infinite product representation {the Weierstrass product 
formula) of the Vigneras multiple gamma function by considering the classical 
limit of the multiple g-gamma functions. In order to get these results, we will use 
the method developed in j22j . Namely, by making use of the Euler-MacLaurin 
summation formula, we derive the Euler-MacLaurin expansion of the multiple 
g-gamma function. Taking the classical limit, we lead to the Euler-MacLaurin 
expansion of the Vigneras multiple gamma function. The higher Stirling formula 
and the Weierstrass product formula follow from this expansion formula 

The details of the proof will be published in the forthcoming paper . 



2 A survey of the multiple gamma function and 
the multiple g-gamma function 

2.1 The gamma function 

The following are well-known facts in the classical analysis: The Bohr-Morellup 
theorem says that the gamma function T{z) is characterized by the three con- 
ditions, 

1. r(z+ 1) = zT{z), 

2. r(i) = i, 

3. ^logr(z + 1) > for z > 0. 

The gamma function is meromorphic on C, and has an infinite product repre- 
sentation 



r(z + i) = e--n 1 + - - ' (1) 



n 

n=l ^ 

where 7 is the Euler Constant. This is usually called the Weierstrass product 
formula. 

The gamma function has an asymptotic expansion, which is called the Stir- 
ling formula, 

logr(z + 1)^ (z + l) log(z -t- 1) - (z + 1) - c'(o) 



E 



2 

B2r 1 



^[2r]2 {z + l)^r-i' 
as z ^ 00 in the sector Ag := {z G C|| argz| < tt — 6} {0 < S < n), where 

ze*- _^Bn{t) 



-z 



n=0 
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Bk = Bk{0) (the Bernoulli number), C(s) is the Riemann zeta- function, C{s) = 
£C{s) and [x]r ^ x{x - 1) ■ ■ ■ {x - r + I). Note that C'(0) = -logv^. 

2.2 The Barnes G-function 

Barnes introduced the function G{z) which satisfies 

1. G{z + 1) =T{z)G{z), 

2. G(l) = l, 

3. ^logG(z + l)>Oforz>0, 

and he called this "G-function" . He proved that the G-function has an infinite 
product representation. 

Giz 4- 1) = e-^C'(0)-4.-^ n { (l + D ' + S) } ■ 

and an asymptotic expansion 

logG{z + l) ^ (^^ - ^yog{z + l)-lz^-^ + ^ + zCiO)~\ogA + Oi^) (2) 

as 2; oo in the sector where A is called the Kinkclin constant. Voros 
showed this constant can be written with the first derivative of the Riemann 
zeta-function (cf ||) 

logA = -C'(-l) + ^. 

2.3 The Vigneras multiple gamma function 

As a generalization of the gamma function and the G-function, Vigneras 
introduced a hierarchy of functions which satisfy 

1. G„(z + 1) = G„_i(z)G„(z), 

2. G„(l) = l, 

3. ^^logG„(z + l) > for z>0, 

4. Go(z)-z. 

and she called these functions "the multiple gamma functions". Applying 
Dufresnoy and Pisot's results P|, she showed that these functions satisfying the 
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above properties are unique and that it has an infinite product representation: 



Gn{z + 1) = exp 



X n 

mGN"-ixN* . 



n-1 



(-1)" (n-l 



s(m) 



exp 



E 



i-i 



n — I \s(m) 



z 



(3) 



where 

En{z) 



E 

m6N'*-i xN* 



E 

k (=0 



(-1)"" 

n — / \^ s(m) 



n — l 



+ (-l)"log 1 + 



s(m) 



:=^logG„-i(z + l), 
:= l'' + 2" + ••• + (z-l)^ 



s(m) := TOi + TO2 
and N* N- {0}. 



m„ for m = (mi,m2,- • -TOn), 



2.4 The g-gamma function 

Throughout this article, we suppose Q < q < 1. A g-analogue of the gamma 
function was defined by Jackson [ p^ : 

1 - 



r(z + l;q) = (l-g)--n(- 



fc=l 

Askey Q pointed out that this function satisfies a g-analogues of the Bohr- 
MoreUup theorem. Namely, V{z\q) satisfies 

1. r(z+l;g) = [z]r(z;g), 

2. r(l;g) = l, 

3. ^logr(z + > for z > 0, 

where [z] := (1 — — q)- 

As q tends to 1 — 0, r(z; q) converges T{z) uniformly with respect to z. A 
rigorous proof of this fact was given by Koornwinder 



mly 
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Inspired by Moak's works [Q, the authors |£2| derived a representation of 
the g-gamma function 

logr(z:g) - (^^l)log(^-:^]+\ogq [\-^di 



2fe-l 



where 



R2m{z;q) (4) 



and the polynomial Pn{x) is defined by the recurrence formula 

Pi(z) = 1, (a;^ - x)-^Pn{x) + nxPnix) = F„+i(a;). 

Each term of the formula converges uniformly as q ^ 1 — 0. So we get 
another proof of the uniformity of the classical limit of logr(z; q). 

2.5 The multiple g-gamma function 

Recently, one of the authors construct the function G„(z; q) which satisfies 
a g-analogue of the generalized Bohr-Morellup theorem 

1. G„(z + 1; g) = G„-i(z; g)G„(z; g), 

2. G„(l;(z) = l, 

3. ^logG„+i(z + l;g) >Oforz>0, 

4. Go(z;g) = [z]. 

We call it "the multiple g-gamma function" . It is given by the following infinite 
product representation: 

Go(2; + l;g) := [z + 1] 



k ( ) 

G„(. + l;g) :=(!-,)-» nU^Y-S") (1 " '^')^"^^''^ [ (5) 



k=l 
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where 



gn{z;u) = 

















\n~\) 



In the next section, we derive a representation of the muhiple g-gamma 
function like (Q) and consider its classical limit. This limit formula gives some 
important properties of the multiple gamma function. 



3 The Euler-MacLaurin expansion 

Our main tool is the Euler-MacLaurin summation formula 

Y.f(')= + ^|f{/('=-i)(iv)^/('^-i)(M)} 



r=M 



k=l 



+ (-1)" 



M 



where Bn{t) = Bn{t — [t]) and [t] denotes an integral part of t. We apply this 
to logG„(z + \;q). 

From the infinite product representation (|^), we have 



\ogGn{z + 



k=l 



-k 



iog(i - g) - E ( „ _ 1 ) " 9^^') 



Tl-l 



+ EG„,,(z){5]Flog(l-<z'=)}, 



where 



3=0 

z — u 
n 



k=l 



)n — l 



By the Euler-MacLaurin summation formula, we obtain 



(6) 



E 

fc=i 



-k 

n — 1 



log(l - 9^+'=) 



(it 



r-l 



n — 1 



t=i 



E — r^n,'--l(2; <?) + I^n,rn{z; q), 



(7) 
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where 



Lr{z) := f , Li(z) :--log(l-g^), L,(l), 
log q 



OO J, 



Lir{z) :— — (Eulcr's polylogarithm) , 

1 - q^+* 



fc=i 



R„,m{z;q) := 











_ (- 


_ j^'^rn— 1 




m! 



n - 1 



log 



1-q 
d" 



z+l 



t=l 



log 



1 - 

1 - 



and „5j is the Stirling number of the first kind defined by 

n 

nSjU^ — u{u — l)(u — 2) • • • (u — n + 1). 

Similarly, we have 

fc=l r=0 r=l 

where 



n+l 



— - E TT-^J+i^'-i^^) 



Bn+i{i) 



P log 



1-9* 
1-g 



dV 



—r < log 



1 -g* 



Substituting (0) and (|) to (|), we have 

logG'„(z + l;q) 

= (the terms containing Lr{z) and i^) 

+ (the term which converges as q ^ 1). 

Lr{z) and Lj. give rise to a divergent part in the expression (^) as g ^ 1 
But, thanks to the identity 



Li+i{z) = — log 



1 -g^ 

1-g 
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■E 



{i-ry. 



^ {l-r)\ r! 



we can show that the first term of (^) vanishes. Therefore, log G„ (z + 1 ; g) reads 
as follows: 

Proposition 1 Suppose SRz > — 1 and m > n. Then we have 
logG„(z + l]q) 



z + l 
n 



Br ( d 



! V dz 



»-0 -(W 



d \ I z 



dz 



n — 1 



■E' 

n— 1 m ^ 

G„j(z)Cj(g) + ^ — [^i^„,r_i(z;g) - Rn^,n{z;q)- 

j=0 r=l 



4 Classical limit 



Now, let us calculate the classical limit. As g ^ 1 — 0, we see that for 5Rz > — 1, 

'1 -q'+^' 



(i) lo, 
(ii) 



1-g 



log (z + l), 



-r-i rc^C^-/ — j-rf^ = — j-{(2 + l) -1}, 



(iii) Fn^r^i{z;q) Fn,r~l{z), 



(iv) i? (z; q) ^ R 

n.rn 



(v) Q(g)->Q:=-5:-f - {tnogO 



(-1)" 

(n + 1)! 



t=i 

n+l 



S„+i(<)f-j {t^logijdt. 



where 



dV-^ \ / -t 



z + t 
n— 1/ + l 
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m! J I J IV'^~1/ V-z + l 

Furthermore, from Vitali's convergence theorem, we can show that the conver- 
gence is uniform on any compact set in {^z > —1}. 

The constant Cj is relevant to the Riemann zeta-function. Indeed, we obtain 

C,=eM-C'i-j))-j-^ (10) 

by applying the Euler-MacLaurin summation formula. 

Thus we have proved the existence of the classical limit of Gn{z + l',q), and 
put it to be 

Gn{z + 1):= \im Gn{z + l;q). 

g— >1— 

The uniformity of the convergence ensures that Gn{z + 1) satisfies 
1. G„(2 + l) = G„_i(0)G„(2), 



2. l^—j logG„(z + l)>0 for z>0, 

3. G„(l) = l, 

4. Go(2 + l) = ^ + l. 

Prom the uniqueness of the function satisfying these conditions, we have 

G„(^-|-l) = G„(^; + 1) for ^z > -1. 
Namely, as q ^ 1 — 0, 

Gn{z + l;q)^Gn{z + l) in {3?^ > -1}. 

Using the functional equation Gn{z + 1) = G„_i(2:)G„(2;), we can show the 
following theorem. 

Theorem 2 As q ^ 1 — 0, Gn{z + 1; (?) converges Gn{z + 1) uniformly on any 
compact set in the domain C\Z<o, and 



\ogGn{z + l) 



z + l\ v-^ Br I a \ I z 



r=l 
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n-1 , 1 1 °° R 



J=0 
^7i.m (^) ■ 

5 The higher StirUng formula 

As |z| — > oo in A5, we can see that 

\Rn,rn{z)\ = 0(z"-'"+l), |F„,,,_i(z)| - 0(z-'-+"). 

Thus, we have proved the higher Stirhng formula: 
Theorem 3 As \z\ 00 in the sector A^, we obtain 
logG'„(z + 1) 



r! V dz / \n — 1 

r— 1 ^ ' ^ 



^ \ '"-1 / \ 1 1 

a \ z \ \ 1 



n-l J- s 00 „ 

Examples of the higher Stirling formula. Let us show some examples. 
In the case that n — I, this formula coincides with the Stirling formula. In each 
case that n = 2, it coincides with (||). Namely, we have 

logG2(z + l) 

- 1) log(z + 1) - - f + i - .c'(0) + C'(-i) 



12z + l ^ [2r]3 (2 + ^ ^ 

In the case that n — S and n = 4, we obtain 
logG3(z + l) 

'Z^ 1\ , , X 11 3 5 n z 13 

\ log(z + 1) z-^ H H 

6 4 24 / ^ 36 24 3 72 
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c'(o) + ^c'(-i)-k'(-2) 



1 1 

r=2 



+ — + ^{z^ ~ (6r - n)z + (4r2 - 16r +16)} 



logG4(z + l) 

log(z + 1) 



' V24 ^ y ~ 720 

44 23 z2 11 31 
'72^ + 9^' 36^+ 144 



12z+l 720 (z + 1)3 \^ 2 2 

+ E rf T r a, ^2.^1 - (12^ - 27)^' + (20r2 - 94r + 111). 



r=3 



-(Sr^ - 56r2 + 134r - 109)} . 

6 The Weierstrass product representation of the 
multiple gamma function 

In this section, we derive the infinite product representation of the muitiple 
gamma function more explicitly than (^. First, we observe the following propo- 
sition. 



Proposition 4 

oc 

exp(C'(-j)) = exp(F,(l)) H U ^ + T ) ^ ' " exp (P,(fc + 1) - P,{k)) 



k=l 

where 

i+i 



Br.(k + 1) 
1 \ ^ + 1 



r=0 
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and the infinite product converges absolutely. 

This proposition is proved by using the Euler-MacLaurin summation formula 
and @. 

Similar calculation shows 

Proposition 5 

n-l 

logG„(^ + 1) = E GnA^)KAz) (11) 

where 

K,{z) ^^'+^^,\+^^ log(z + 1) - C'(-j) + + 1) 



J + 1 

OC p 



k=l 



B.j+i{z + k + l) , fz + k + l 
Pj{z + k + l)-Pj{z + k) + '- log ' 



j + 1 °\ z + k 



Furthermore, the infinite sum of the last term is absolute convergent. 

By using Proposition 4 and careful consideration on Kj{z), we can see 

■J 1 ..-l-r/-ll \ ^ 



E + fc)'^ log (i + f ) + E f:) ^^-^ E ^^^fc"^-' 



where 



k=l (. r=0 ^ ^ 1 = 1 



J 

J 



r=0 



The infinite sum in this formula converges absolutely. 
Finally we can derive the following theorem: 
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Theorem 6 For n gN, we have 



Gn{z + 1) = exp n (l + I) exp ($„(^, k)) I , 



where 



j=0 r=0 



1 (^AV ('^-'' 

r\ \du J \n — 1 



X C'(-r) 



n- 1 



Examples of the Weierstrass product representation. By making 

use of Theorem 6, the Weierstrass product representation of the multiple gamma 
function is derived explicitly. We give some examples. In the case that n = 1 
and n = 2, the results coincide with the Weierstrass representation of the gamma 
function and the G-function. In the case that n = 3 and n = 4, we obtain 



{ z"^ 7 , ziz - 1) , , ~ 



n 

fe=i 



fc(fc+i) 



-k 



Gi{z + 1) 



61 4 13 o 19 2 5 
^"P^Tii" +18^ +144^ -24^ 

_£c'(-2) + ^C'(-i) - ^^^^^C'(0) - £l^i£!±i£!7 



n 

fe=i 



fc(fc+l)(fc + 2) 



exp 



1 



24 6 6 / fc 



18 



:_ _ r ) - lfc2 
12 2 7 6 
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